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Abstract 

The NA48/2 experiment at the CERN SPS carried out data taking in 2003 
and 2004. Analysis of the selected data samples of 7,146 — + 7r ± e + e _ 
decay candidates with 0.6% background, 1,164 K — > 7r ± 77 candidates with 
3.3% background, and 120 .ff — > 7r ± 7e + e~ candidates with 6.1% background 
allowed precise measurements of branching fractions and other characteristics 
of these rare kaon decays. 



Introduction 



Radiative nonleptonic kaon decays represent a source of information on the 
structure of the weak interactions at low energies, and provide crucial tests 
of the Chiral Perturbation Theory (ChPT). The current paper presents new 
results related to study of the — > 7r ± e + e~, — > 71^77, and K — > 
Tr^e+e" decays by the NA48/2 experiment at the CERN SPS. 

The flavour-changing neutral current process — ► 7r ± e + e _ , induced 
at one-loop level in the Standard Model and highly suppressed by the GIM 
mechanism, has been described by the ChPT 03 ; several models predicting the 
form factor characterizing the dilepton invariant mass spectrum and the decay 
rate have been proposed ^ © The decay is fairly well explored experimen- 
tally: it was first studied at CERN © , followed by BNL E777® and E865© 
measurements . 

The — > 7r ± 77 and — ► 7r ± 7e + e _ decays similarly arise at one-loop 
level in the ChPT. The decay rates and spectra have been computed at leading 
and next-to- leading orders^ and strongly depend on a single theoretically 
unknown parameter c. The experimental knowledge of these processes is rather 
poor: before the NA48/2 experiment, only a single observation of 31 K — > 

7r ± 77 candidates was made^, while the — > n :iz je + e~ decay was not 
observed at all. 

The paper is organized as follows. In Section 1, a description of the 
NA48/2 experiment is given. Section 2 is devoted to a rather detailed de- 
scription of the — ► 7r ± e + e~ analysis and its preliminary results, which is 
the main topic of the paper. Section 3 briefly presents the preliminary results 
of the — > 7r ± 77 analysis; a more detailed discussion is reserved for the 
Moriond QCD 2008 conference. Section 4 briefly presents the final results of 

the — > 7r ± 7e + e~ analysis, which have recently been published 0^ . Finally 
the conclusions follow. 



1 The NA48/2 experiment 

The NA48/2 experiment, designed to excel in charge asymmetry measure- 
ments j i s based on simultaneous K + and K~ beams produced by 400 
GeV/c primary SPS protons impinging at zero incidence angle on a beryllium 
target of 40 cm length and 2 mm diameter. Charged particles with momentum 
(60 ± 3) GeV/c are selected by an achromatic system of four dipole magnets 
with zero total deflection ('achromat'), which splits the two beams in the ver- 
tical plane and then recombines them on a common axis. Then the beams 
pass through a defining collimator and a series of four quadrupoles designed 
to produce focusing of the beams towards the detector. Finally the two beams 




Figure 1: Schematic lateral view of the NA48/2 beam line (TAX17,18: motor- 
ized beam dump /collimators used to select the momentum of the K + and K~ 
beams; FDFD/DFDF: focusing set of quadrupoles, KABES1-3: kaon beam 
spectrometer stations), decay volume and detector (DCH1-4-' drift chambers, 
HOD: hodoscope, LKr: EM calorimeter, HAC: hadron calorimeter, MUV: 
muon veto). The vertical scales are different in the two parts of the figure. 



are again split in the vertical plane and recombined in a second achromat. The 
layout of the beams and detectors is shown schematically in Fig. [T] 

The beams then enter the decay volume housed in a 114 m long cylindrical 
vacuum tank with a diameter of 1.92 m for the first 65 m, and 2.4 m for the 
rest. Both beams follow the same path in the decay volume: their axes coincide 
within 1 mm, while the transverse size of the beams is about 1 cm. With 
7 x 10 11 protons incident on the target per SPS spill of 4.8 s duration, the 
positive (negative) beam flux at the entrance of the decay volume is 3.8 x 10 7 
(2.6 x 10 7 ) particles per pulse, of which 5.7% (4.9%) are K+ (R-). The 
K + 1 K~ flux ratio is about 1.8. The fraction of beam kaons decaying in the 
decay volume at nominal momentum is 22%. 

The decay volume is followed by a magnetic spectrometer housed in a tank 
filled with helium at nearly atmospheric pressure, separated from the vacuum 
tank by a thin (0.31%Ao) Kevlar composite window. A thin-walled aluminium 
beam pipe of 16 cm outer diameter traversing the centre of the spectrometer 
(and all the following detectors) allows the undecayed beam particles and the 
muon halo from decays of beam pions to continue their path in vacuum. The 
spectrometer consists of four drift chambers (DCH): DCH1, DCH2 located 
upstream, and DCH3, DCH4 downstream of a dipole magnet. The magnet 
provides a horizontal transverse momentum kick Ap — 120 MeV/c for charged 
particles. The DCHs have the shape of a regular octagon with a transverse 



size of about 2.8 m and a fiducial area of about 4.5 m 2 . Each chamber is 
composed of eight planes of sense wires arranged in four pairs of staggered 
planes oriented horizontally, vertically, and along each of the two orthogonal 
45° directions. The spatial resolution of each DCH is a x = <j y = 90 /im. The 
nominal spectrometer momentum resolution is o v jp = (1.02 © 0.044 ■ p)% (p in 
GeV/c). 

The magnetic spectrometer is followed by a plastic scintillator hodoscope 
(HOD) used to produce fast trigger signals and to provide precise time mea- 
surements of charged particles. The hodoscope has a regular octagonal shape 
with a transverse size of about 2.4 m. It consists of a plane of horizontal and 
a plane of vertical strip-shaped counters. Each plane consists of 64 counters 
arranged in four quadrants. Counter widths (lengths) vary from 6.5 cm (121 
cm) for central counters to 9.9 cm (60 cm) for peripheral ones. 

The HOD is followed by a liquid krypton electromagnetic calorimeter 

(LKr) use d for photon detection and particle identification. It is an al- 
most homogeneous ionization chamber with an active volume of 7 m 3 of liquid 
krypton, segmented transversally into 13248 projective cells, 2x2 cm 2 each, 
by a system of Cu— Be ribbon electrodes, and with no longitudinal segmenta- 
tion. The calorimeter is 27Xq deep and has an energy resolution a(E)/E = 
0.032/\/Ee0.09/£;e0.0042 (E in GeV). Spatial resolution for a single electro- 
magnetic shower is a x — a y — 0.42/-\/E©0.06 cm for the transverse coordinates 
x and y. 

The LKr is followed by a hadronic calorimeter (HAC) and a muon detector 
(MUV), both not used in the present analysis. A detailed description of the 

components of the NA48 detector can be found elsewhere . The NA48 / 2 
experiment took data during two runs in 2003 and 2004, with about 60 days 
of effective running each. About 18 x 10 9 events were recorded in total. 

In order to simulate the detector response, a detailed GEANT-based H^JI 
Monte Carlo (MC) simulation is employed, which includes full detector geome- 
try and material description, stray magnetic helds, DCH local inefficiencies and 
misalignment, detailed simulation of the kaon beam line, and time variations 
of the above throughout the running period. Radiative corrections are applied 
to kaon decays using the PHOTOS package 

2 K^ 1 — > Tr^e+e - analysis 

The — > 7r ± e + e~ rate is measured relatively to the abundant K — > n^n^ 
normalization channel (with tt^ — > e + e~-f). The final states of the signal and 
normalization channels contain identical sets of charged particles. Thus elec- 
tron and pion identification efficiencies, potentially representing a significant 
source of systematic uncertainties, cancel in the first order. 



2.1 Event selection 



Three- track vertices (compatible with the topology of if — > n e + e and 

— > Tr^irp decays) are reconstructed using the Kalman filter algorithm US' 
by extrapolation of track segments from the upstream part of the spectrometer 
back into the decay volume, taking into account the measured Earth's mag- 
netic field, stray field due to magnetization of the vacuum tank, and multiple 
scattering in the Kevlar window. 

A large part of the selection is common to the signal and normalization 
modes. It requires a presence of a vertex satisfying the following criteria. 

• Total charge of the three tracks: Q = ±1. 

• Vertex longitudinal position is inside fiducial decay volume: Vertex > 

■^final collimator- 

• Particle identification is performed using the ratio E/poi track energy de- 
position in the LKr to its momentum measured by the spectrometer. The 
vertex is required to be composed of one pion candidate (E/p < 0.85), 
and two opposite charge e^- candidates (E/p > 0.95). No discrimination 
of pions against muons is performed. 

• The vertex tracks are required to be consistent in time (within a 10 ns 
time window) and consistent with the trigger time, to be in DCH, LKr and 
HOD geometric acceptance, and to have momenta in the range 5 GeV/c < 
p < 50 GeV/c. Track separations are required to exceed 2 cm in the 
DCH1 plane to suppress photon conversions, and to exceed 15 cm in the 
LKr plane to minimize particle misidentification due to shower overlaps. 

If multiple vertices satisfying the above conditions are found, the one with the 
best fit quality is considered. The following criteria are then applied to the 
reconstructed kinematic variables to select the K — > 7r ± e + e~ candidates. 

• 7r ± e + e _ momentum within the beam nominal range: 54GeV/c < IpVeel < 
66 GeV/c. 

• ir ± e + e~ transverse momentum with respect to the measured beam tra- 
jectory: <p% < 0.5 x 10~ 3 (GeV/c) 2 . 

• n ± e+e'~ invariant mass: 475 MeV/c 2 < M xee < 505 MeV/c 2 . 

• Suppression of the K — » w background defining the visible kine- 
matic region: z — (M ee /Mic) 2 > 0.08, which approximately corresponds 
to M ee > 140 MeV/c 2 . 




Figure 2: Left: reconstructed spectrum o/7r ± e + e _ invariant mass; data (dots) 
and MC simulation (filled area). Right: the computed dY^ ee /dz (background 
subtracted, trigger efficiencies corrected for) and the results of fits according to 
the considered models. 



Independently, a presence of a LKr energy deposition cluster (photon candi- 
date) satisfying the following principal criteria is required to select the i"T — > 
7r ± 7r^ ) candidates. 

• Cluster energy E > 3 GeV, cluster time consistent with the vertex time, 
sufficient transverse separations from track impact points at the LKr 
plane (R^ > 30 cm, i? e7 > 10 cm). 

• e + e~7 invariant mass compatible with a ir° decay: \M eei — M„o\ < 
10 MeV/c 2 . 

• The same conditions on reconstructed 7r ± e + e~7 total and transverse mo- 
menta as used for ir ziz e + e~ momentum in the — ► 7r ± e + e _ selection. 

• n ± e + e~-f invariant mass: 475 MeV/c 2 < M 7ree7 < 510 MeV/c 2 . 
2.2 Signal and normalization samples 

The reconstructed ir ziz e + e~ invariant mass spectrum is presented in Fig. [5] (left 
plot). The 7r ± e + e _ mass resolution is ov ee = 4.2 MeV/c 2 , in agreement with 
MC simulation. The e + e~ mass resolution computed by MC simulation is 
a ee = 2.3 MeV/c 2 . 



In total 7,146 K^- — ► 7r ± e + e _ candidates are found in the signal region. 
After the kinematical suppression of the ifi D decays, residual background con- 
tamination mostly results from particle misidcntification (i.e. e ± identified as 
71"^ and vice versa). The following relevant background sources were identi- 
fied with MC simulations: (1) — > ^it^ with misidentified and tt^; 
(2) — > 7r^ ) e ± ^ with a misidentified e ± from the 77^ decay. Background 

estimation by selecting the strongly suppressed 03) lepton number violating 
— > 7r T e ± e ± ("same-sign") candidates was considered the most reliable 
method. For the above two background sources, the expected mean numbers 
and kinematic distributions of the selected same-sign candidates are identical to 
those of background events (up to a negligible acceptance correction) . In total 
44 events pass the same-sign selection, which leads to background estimation 
of (0.6 ± 0.1)%. This result was independently confirmed with MC simulation 
of the two background modes. 

In total 12.228 x 10 6 — > n^n^ candidates are found in the signal re- 
gion. The only significant background source is the semileptonic — > tt° d ii ± v 
decay. Its contribution is not suppressed by particle identification cuts, since 
no 7r//i separation is performed. The background contamination is estimated 
to be 0.15% by MC simulation. 

2.3 Trigger chain and its efficiency 

Both — ► 7r ± e + e _ and — ► 7r ± 7r^ ) samples (as well as K ± — ► "Stt^) 
are recorded via the same two-level trigger chain. At the first level (LI), a 
coincidence of hits in the two planes of the HOD in at least two of the 16 non- 
overlapping segments is required. The second level (L2) is based on a hardware 
system computing coordinates of hits from DCH drift times, and a farm of 
asynchronous processors performing fast track reconstruction and running a 
selection algorithm, which basically requires at least two tracks to originate 
in the decay volume with the closest distance of approach of less than 5 cm. 
LI triggers not satisfying this condition are examined further and accepted 
nevertheless if there is a reconstructed track not kinematically compatible with 
a tt^tt decay of a having momentum of 60 GeV/c directed along the beam 
axis. 

The NA48/2 analysis strategy for non-rare decay modes involves direct 
measurement of the trigger efficiencies using control data samples of downscaled 
low bias triggers collected simultaneously with the main triggers. However di- 
rect measurements are not possible for the — > TT^e + e~ events due to very 
limited sizes of the corresponding control samples. Dedicated simulations of 
LI and L2 performance (involving, in particular, the measured time dependen- 
cies of local DCH and HOD inefficiencies) were used instead. The simulated 
efficiencies and their kinematic dependencies were compared against measure- 
ments for the abundant — > ir^ir^ and K — > Tt^it+Ti - decays in order to 



validate the simulations. 

The simulated values of LI and L2 inefficiencies for the selected — ► 
^7^% sample are Eli — 0.37%, Eli — 0.80%. The values of the integral trigger 
inefficiencies for the K ± — ► ir :il e + e~ sample depend on the a priori unknown 
form factor; the corrections are applied differentially in bins of dilepton in- 
variant mass. Indicative values of inefficiencies computed assuming a realistic 
linear form factor with a slope 8 = 2.3 are Eli = 0.06%, £L2 = 0.42%. The 
— > 7T^~ 7r^j sample is affected by larger inefficiencies due to a smaller invari- 
ant mass of the e + e~ system, which means that the leptons are geometrically 
closer. 

2.4 Theoretical input 

The decay is supposed to proceed through one photon exchange, resulting in a 
spectrum of the z — {M ee /MK) 2 kinematic variable sensitive to the form factor 

W(z)^: 

wherer e = m e /MK, r n = iti^/Mk, and X(a,b,c) = a 2 + b 2 +c 2 — 2ab—2ac—2bc. 
On the other hand, the spectrum of the angle O^e between ir and e + in the e + e~ 
rest frame is proportional to sin 2 6* we , and is not sensitive to W[z). 

The following parameterizations of the form factor are considered 

in the present analysis. 

1. Linear: W(z) — GfM|-/o(1 + Sz) with free normalization and slope 
(fo,S). 

2. Next-to-leading order ChPT©: W(z) = G F M 2 i {a+ + b+z) + W™(z) 
with free parameters (a+, b + ), and an explicitly calculated pion loop term 
W WK {z). 

3. The Dubna version of ChPT parameterization involving meson form fac- 
tors: W(z) = W(M a , M p , z)©, with resonance masses (M a , M p ) treated 
as free parameters. 

The goal of the analysis is extraction of the form factor parameters in the 
framework of each of the above models, and computation of the corresponding 
branching fractions BRi2,3. 



2.5 Fitting procedure 



The values of dY^ee/dz in the centre of each i-bin of z, which can be directly 
confronted to the theoretical predictions ((I}, are then computed as 

(dT^/dz), = ^-^l ■ A ^' £2 ^ ■ BR(^ ± - vr^ ) • BR(^) ■ (2) 
N 27I Ai{l-£i) Az 

Here JVi and Nf are the numbers of observed — > n ± e + e^ candidates and 
background events in the i-th bin, N 2w is the number of — ► 7^7^ events 
(background subtracted), A4 and e< are geometrical acceptance and trigger 
inefficiency in the i-th bin for the signal sample (computed by MC simulation) , 
A 2 tt = 2.94% and e 2v = 1.17% are those for — > events, Tk is the 

nominal kaon width US', Az is the chosen width of the z bin, BR(_ft' ± — > 

t^tt ) = (20.64±0.08)% (FlaviaNet average®), BR(w° D ) = (1.198 ± 0.032)% 

(PDG average®). 

The computed values of dY vee jdz vs z are presented in Fig. [2] (right plot) 
along with the results of the fits to the three considered models. BR(if ± — > 
7r ± e + e _ ) in the full kinematic range corresponding to each model are then 
computed using the measured parameters, their statistical uncertainties, and 
correlation matrices. 

In addition, a model-independent branching fraction BR m ; in the visible 
kinematic region z > 0.08 is computed by integration of dT^ee/dz. BR m ; is to 
a good approximation equal to each of the model-dependent BRs computed in 
the restricted kinematic range z > 0.08. 

2.6 Systematic uncertainties 

The following sources of systematic uncertainties were studied. 

1. Particle identification. Imperfect MC description of electron and pion 
identification inefficiencies f e and f w can bias the result only due to the mo- 
mentum dependence of the inefficiencies, due to identical charged particle com- 
position, but differing momentum distributions of the signal and normalization 
final states. Inefficiencies were measured for the data to vary depending on 
particle momentum in the ranges 1.6% < f n < 1.7% and 1.1% < f e < 1.7% 
in the analysis track momentum range. Systematic uncertainties due to these 
momentum dependencies not perfectly described by MC were conservatively 
estimated assuming that MC predicts momentum- independent f e and f % . 

2. Beam line description. Despite the careful simulation of the beamlinc 
including time variations of its parameters, the residual discrepancies of data 
and MC beam geometries and spectra bias the results. To evaluate the related 
systematic uncertainties, variations of the results with respect to variations of 
cuts on track momenta, LKr cluster energies, total and transverse momenta of 



Parameter 


e, 7r 


Beam 


Background 


Trigger 


Rad. 


Fitting 




ID 


spectra 


subtraction 


efficiency 


corr. 


method 


S 


0.01 


0.04 


-0.04 ± 0.04 


-0.03 ± 0.03 


0.05 


0.03 


fa 


0.001 


0.006 


0.002 ± 0.002 


0.000 ± 0.001 


0.006 


0.003 


a+ 


0.001 


0.005 


-0.001 ± 0.001 


-0.001 ± 0.002 


0.005 


0.004 


b+ 


0.009 


0.015 


0.017 ± 0.017 


0.016 ± 0.015 


0.015 


0.010 


M a /GeV 


0.004 


0.009 


0.008 ± 0.008 


0.006 ± 0.006 


0.009 


0.006 


Af fc /GeV 


0.002 


0.003 


0.003 ± 0.003 


0.003 ± 0.003 


0.004 


0.002 


BRi,2, 3 xlO Y 


0.02 


0.02 


-0.01 ± 0.01 


-0.02 ± 0.01 


0.01 


0.02 


BR mi xl0 7 


0.02 


0.01 


-0.01 ± 0.01 


-0.02 ± 0.01 


0.01 


n/a 



Table 1: Summary of corrections and systematic uncertainties (excluding the 
external ones). 



the final states ir ± e + e~('j), and track distances from beam axis in DCH planes 
were studied. 

3. Background subtraction. As discussed above, the same-sign event spec- 
trum is used for background estimation in the Tr^e + e^ sample. The method 
has a limited statistical precision (with an average of 2 same-sign event in a bin 
of z). Furthermore, the presence of the component with two e + e~ pairs (due 
to both ir^ decays and external conversions) with a non-unity expected ratio 
of same-sign to background events biases the method. The uncertainties of the 
measured parameters due to background subtraction were conservatively taken 
to be equal to the corrections themselves. 

4. Trigger efficiency. As discussed earlier, the corrections for trigger 
inefficiencies were evaluated by simulations. In terms of decay rates, LI and 
L2 corrections have similar integral magnitudes of a few 10 -3 . No uncertainty 
was ascribed to the LI correction, due to relative simplicity of the trigger 
condition. On the other hand, the uncertainty of the L2 efficiency correction 
was conservatively taken to be equal to the correction itself. 

5. Radiative corrections. Uncertainties due to the radiative corrections 
were evaluated by variation of the lower 7r ± e + e _ invariant mass cut. 

6. Fitting method. Uncertainties due to the fitting procedure were eval- 
uated by variation of the z bin width. 

7. External input. Substantial uncertainties arise from the external input, 
as BR(7r ± 7r^ ) ) is experimentally known only with 2.7% relative precision , 
The only parameter not affected by an external uncertainty is the linear form 
factor slope S describing only the shape of the spectrum. 

The applied corrections and the systematic uncertainties (excluding the 
external ones presented later) are summarized in Table [TJ 



5 = 

fo = 
BRi x 10 7 = 


2.35 ± 0.15 sta t. 
0.532 ± 0.012 stat . 
3.02 ± 0.04 stat . 


± 0.09 syst . 
± 0.008 syst . 
± 0.04 syst . 


± o.oo cxt . 

± 0.007 oxt . 
± 0.08 cxt . 


= 2.35 ± 0.18 
= 0.532 ± 0.016 
= 3.02 ± 0.10 


a+ = 

b+ = 
BR 2 x 10 7 = 


-0.579 ± 0.012 stat . 
-0.798 ± 0.053 sta t. 
3.11 ± 0.04 stat . 


± 0.008 syst . 
± 0.037 syst . 
± 0.04 syst . 


± 0.007 oxt . 
± 0.017 cxt . 
± 0.08 cxt . 


= -0.579 ± 0.016 
= -0.798 ± 0.067 
3.11 ± 0.10 


M a /GeV = 
M p /GeV = 
BR 3 x 10 7 = 


0.965 ± 0.028 sta t. 

0.711 ± 0.010 rt at. 

3.15 ± 0.04 stat . 


± 0.018 syst . 
± 0.007 syst . 
± 0.04 syst . 


± 0.002 oxt . 
± 0.002 ext . 
± 0.08 cxt . 


= 0.965 ± 0.033 
= 0.711 ± 0.013 
= 3.15 ± 0.10 


BR mi x 10 Y = 


2.26 ± 0.03 stat . 


± 0.03 syst . 


± 0.06 cxt . 


2.26 ± 0.08 



Table 2: Results of fits to the three considered models, and the model- 
independent BR mi (z > 0.08). 



2.7 Results and discussion 

The measured parameters of the considered models and the corresponding BRs 
in the full z range, as well the model-independent BR m i(z > 0.08), with their 
statistical, systematic, and external uncertainties are presented in Tabled The 
correlation coefficients between the pairs of model parameters, not listed in the 
table, are p(S, f ) = -0.963, p(a+, b+) = -0.913, and p(M a ,M p ) = 0.998. 

Fits to all the three models are of reasonable quality, however the linear 
form-factor model leads to the smallest \ 2 ■ The data sample is insufficient to 
distinguish between the models considered. 

The obtained form factor slope S is in agreement with the previous mea- 
surements based on K + — ► ir + e + e~ El ®) and — > tt^ p, + pT samples, 
and, further confirms the contradiction of the data to meson dominance mod- 
els^*. The obtained fo, a+ and fe+ are in agreement with the only previous 
measurement© . The measured parameters M a and M p are a few % away from 

the nominal masses of the resonances . 

The branching ratio in the full kinematic range, which is computed as the 
average between the two extremes corresponding to the models (1) and (3), 
and includes an uncertainty due to extrapolation into the inaccessible region 
z < 0.08, is 

BR=(3.08±0.04 stat .±0.04 syst ±0.08e X t±0.07 modcl )xlO- 7 = (3.08±0.12)xlO- 7 . 

It should be stressed that a large fraction of the uncertainty of this result is 
correlated with the earlier measurements. A comparison to the precise BNL 

E865 measurement ® dismissing correlated uncertainties due to external BRs 
and model dependence, and using the same external input, shows a 1.4(7 dif- 
ference. In conclusion, the obtained BR is in agreement with the previous 



measurements . 

Finally, a first measurement of the direct CP violating asymmetry of K + 
and K~ decay rates in the full kinematic range was obtained by performing 
BR measurements separately for K + and K~ and neglecting the correlated 
uncertainties: A(K± ee ) = (BR+ - BR~)/(BR + + BR") = (-2.1 ± 1.5 stat . ± 
0.3 sys t.)%. The result is compatible to no CP violation. However its precision 
is far from the theoretical expectation© of lA^^,)! ~ 10~ 5 . 

3 K^ 1 — > tv^jj analysis 

The K ± — -> 71^77 rate is measured relatively to the — > ^tt normalization 
channel. The signal and normalization channels have identical particle com- 
position of the final states, and the only cut differing for the two channels is 
the one on the 77 invariant mass. The used trigger chain involves the so called 
"neutral trigger" based on requirement of minimal number of energy deposition 
clusters in the LKr calorimeter. 

About 40% of the total NA48/2 data sample have been analyzed, and 
1,164 K ± — > 7r ± 77 decay candidates (with background contamination esti- 
mated by MC to be 3.3%) are found, which has to be compared with the 

only previous measurement © involving 31 decay candidates. The recon- 
structed spectrum of 77 invariant mass in the accessible kinematic region 
M 77 > 0.2 GeV/c 2 is presented in Fig. [3l along with a MC expectation as- 
suming ChPT 0(p 6 ) distribution © with a realistic parameter c = 2. ChPT 
predicts an enhancement of the decay rate (cusp-like behaviour) at the tttt mass 
threshold m 77 ss 280 MeV/c 2 , independently of the value of the c parameter. 
The observed spectrum provides the first clean experimental evidence for this 
phenomenon. 

As the first step of the analysis, the partial width of the decay was mea- 
sured assuming the ChPT 0(p 6 ) shape with a fixed parameter c = 2. The 
following preliminary result, which is in agreement with the ChPT computa- 
tion for c = 2, was obtained: 

BR = (1.07± 0.04 8tat . ± 0.08 syst .) x lO" 6 . 

A combined fit of the m 77 spectrum shape and the decay rate is foreseen to 
measure the c parameter. 

4 K ^ — y 7r = ' = 7e+e - analysis 

The K — > 7r ± 7e + e _ rate is measured relatively to the — > tt^tt^ normal- 
ization channel. The signal and normalization channels have identical particle 
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Figure 3: The reconstructed spectrum 0/77 invariant mass for the — > 71^77 
decay (dots), and its comparison to MC expectation assuming ChPT 0(p 6 ) 
distribution with c = 2 (filled area). 

composition of the final states. The same trigger chain as for the collection of 
K ± — > 7r ± e + e _ is used. 

With the full NA48/2 data sample analyzed, 120 K^- — > ■K ± ^e + e" decay 
candidates (with the background estimated by MC to be 6.1%) are found in the 
accessible kinematic region M jee > 0.26 GeV/c 2 . This is the first observation 
of this decay mode. The reconstructed spectrum of r ye + e~ invariant mass is 
presented in Fig. [5J along with MC expectations for background contributions. 
The spectrum provides another evidence for the rate enhancement at the tttt 
mass threshold. 

The final results of the analysis have recently been published ^® . The 
model-independent partial width in the accessible kinematic region is measured 
to be 

BR(Af 7ee > 0.26 GeV/c 2 ) = (1.19 ± 0.12 stat . ± 0.04 syst .) x lO" 8 . 

The ChPT parameter c assuming 0(p 4 ) distibution © was measured to be 
c= 0.90 ±0.45. 

Conclusions 

A precise study of the — > 7r ± e + e _ decay has been performed. The data 
sample and precision are comparable to world's best ones, the preliminary 
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Figure 4: The reconstructed spectrum of r ye + e invariant mass for the 
7r ± 7e + e~ decay (dots), and MC background expectations (filled areas). 



results are in agreement with the previous measurements, and the first limit 
on CP violating charge asymmetry has been obtained. 

A precise study of the K ± — > 71^77 has been performed. The first clear 
evidence for a rate enhancement at 7T7r mass threshold has been obtained. The 
preliminary measurement of BR agrees with the ChPT prediction. A detailed 
spectrum shape study is foreseen. 

The first observation of the K ± — > 7r ± 7e + e _ decay, and measurement of 
its parameters, including the BR, have been performed. The M 7ee spectrum 
provides an independent evidence for the cusp at the 7T7t mass threshold. 
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